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Spin-orbit coupling characterizes quantum systems such as atoms, nuclei, hypernuclei, quarkonia, 
etc., and is essential for understanding their spectroscopic properties. Depending on the system, the 
effect of spin-orbit coupling on shell structure is large in nuclei, small in quarkonia, perturbative in 
atoms. In the standard non-relativistic reduction of the single-particle Dirac equation, we derive a 
universal rule for the relative magnitude of the spin-orbit effect that applies to very different quantum 
systems, regardless of whether the spin-orbit coupling originates from the strong or electromagnetic 
interaction. It is shown that in nuclei the near equality of the mass of the nucleon and the difference 
between the large repulsive and attractive potentials explains the fact that spin-orbit splittings are 
comparable to the energy spacing between major shells. For a specific ratio between the particle mass 
and the effective potential whose gradient determines the spin-orbit force, we predict the occurrence 
of giant spin-orbit energy splittings that dominate the single-particle excitation spectrum. 

PACS numbers: 21 . 10 .-k, 21 . 60 .Jz, 21 . 80 .-|-a 


I. INTRODUCTION 

Spin-orbit coupling, known since the early days of 
quantum mechanics, is among the most studied effects 
related to the spin of a particle. After nearly a cen¬ 
tury, spin-orbit effects continue to provide a basis for a 
variety of new phenomena in diverse fields of quantum 
physics, such as spintronics [1], topological insulators [2], 
cold atomic gases [3], atomic nuclei far from stability [4], 
etc. In nuclear physics, in particular, the strong cou¬ 
pling between the orbital angular momentum and spin 
of a nucleon accounts for the empirical magic numbers 
and shell gaps. The energy spacings between spin-orbit 
partner states can be as large as the gaps between major 
shells in atomic nuclei, but the spin-orbit splittings can 
also be considerably reduced in short-lived nuclei with 
extreme isospin values W- Starting from the nuclear 
physics case, we investigate the ratio between the ma¬ 
jor energy spacings of levels characterized by principal 
single-particle quantum numbers and the energy splitting 
of spin-orbit partner states in systems of bound fermions. 
This ratio can differ, of course, by orders of magnitude in 
systems in which the spin-orbit coupling originates from 
different underlying interactions, e.g. the strong effective 
nuclear force or the electromagnetic interaction. It is, 
therefore, interesting to try to find a common character¬ 
istic for the spin-orbit splitting in systems as diverse as 
atoms, nuclei, hypernuclei, quarkonia, etc. 


II. SPIN-ORBIT COUPLING IN ATOMIC 
NUCLEI 

In the nuclear relativistic mean-field framework [7H9] 
a nucleus is considered as a system of independent nu¬ 


cleons moving in local self-consistent scalar and vector 
potentials. The single-nucleon dynamics is governed by 
the Dirac equation: 

[d-p + V + P{m + 5)] V’i = Eitpi (1) 

where tpi denotes the Dirac spinor: 


X* 


( 2 ) 


for the *-th nucleon. For simplicity we only con¬ 
sider spherical nuclei and assume time-reversal symme¬ 
try (pairwise occupied states with Kramers degeneracy), 
which ensures that the only non-vanishing components 
of the vector fields are the time-like ones and thus there 
is no net contribution from nucleon currents. The lo¬ 
cal vector V and scalar S potentials are uniquely deter¬ 
mined by the actual nucleon density and scalar density 
of a given nucleus, respectively. In the ground state A 
nucleons occupy the lowest single-nucleon orbitals, de¬ 
termined self-consistently by the iterative solution of the 
Dirac equation ([^. If one expresses the single-nucleon 
energy as = m + Si, where m is the nucleon mass, 
and rewrites the Dirac equation as a system of two equa¬ 
tions for (f)i and Xi> then, noticing that for bound states 
ei « m, 


X» 


I 

2M{r) 


{^■P) (t'r 


( 3 ) 


to order Ei/m, and 


M{r) = m +-{S{r)— V{r)) . 


( 4 ) 


The equation for the upper component (pi of the Dirac 
spinor reduces to the Schrodinger-like form nniiii] 




= EiC, 


( 5 ) 
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for a nucleon with effective mass A4{r) in the potential 
U{r) = V{r) + S{r), and with the spin-orbit potential: 


yLS ^ 


1 1 d 


2M?‘(r) r dr 




( 6 ) 


The spin-orbit coupling plays a crucial role in nuclear 
structure, and its inclusion in the effective single-nucleon 
potential is essential to reproduce the empirical magic 
numbers. The relativistic mean-field framework, in par¬ 
ticular, naturally includes the nucleon spin degree of free¬ 
dom, and the resulting spin-orbit potential emerges au¬ 
tomatically with the empirical strength [T^]. The nu¬ 
clear spin-orbit potential originates from the difference 
between two large fields: the vector potential V (short- 
range repulsion) with typical strength of « 350 MeV, 
and the scalar potential S (medium-range attraction), 
typically of the order of —400 MeV in nucleonic mat¬ 
ter and finite nuclei. In the context of in-medium QCD 
sum rules m, the strong scalar and vector mean fields 
experienced by nucleons can be associated with the lead¬ 
ing density dependence of the chiral (quark) condensate, 
{qq), and the quark density {q^q}- In the mean-field phe¬ 
nomenology the sum of these two fields V+S « —50 MeV 
provides the confining potential that binds the nucleons 
in a nucleus, whereas the large difference V — S ^ 750 
MeV determines the pronounced energy spacings be¬ 
tween spin-orbit partner states in finite nuclei, of the 
order of several MeV [Ml[mEllis]. A puzzling coin¬ 
cidence that we would like to explore is that the largest 
spin-orbit splittings for intruder states are comparable in 
magnitude to the energy gaps between major shells of 
the nuclear potential. 

The aim of this study is to evaluate the typical ratio 
between the energy spacings of levels characterized by 
principal single-particle quantum numbers and the en¬ 
ergy splitting of spin-orbit partner states (fine structure). 
Even though the values of this ratio span orders of magni¬ 
tude for different bound quantum systems, we will show 
that it is basically governed by two quantities that char¬ 
acterize a given system, irrespective whether the binding 
and spin-orbit potentials originate from the strong (nu¬ 
clei) or electromagnetic (atoms) interactions. 

In the nuclear case the interaction is of short range 
and the self-consistent potentials display a spatial distri¬ 
bution that corresponds to the actual single-nucleon den¬ 
sity. The expression for the spin-orbit potential Eq. ®> 
can, therefore, be rewritten in the following form: 


where 


yLS 


F(r) 


2p(r)r 


Fir) ^ _ 

[m- liV{r) - S'(r))]^ ’ 


( 7 ) 

( 8 ) 


and p(r) denotes the self-consistent ground-state density 
of a nucleus with A nucleons. For a typical approxima¬ 
tion of the single-nucleon potential, such as the harmonic 


oscillator, or the more realistic Woods-Saxon potential, 
one can show [I6j : 


2p(r)r Rl 


( 9 ) 


where i?o = tq » 1.2 fm and, together with < 

I ■ s >= 1/2 for j = 1 + 1/2, and < I ■ s >= —1/2(1 + 1) 
for j = Z — 1/2, the energy spacing between spin-orbit 
partner states can be approximated by: 


|A < > 


Rl 


( 10 ) 


Let us consider the ratio between the major energy spac¬ 
ings and the spin-orbit splitting. For the harmonic oscil¬ 
lator potential one finds m 


fiwo = — 

Un 


h -2Uo 


( 11 ) 


where, in our case, the depth of the potential is Uq = 
t/(r = 0) = V(0) -I- 5(0). Therefore, 


hujo 


X = 


A < > 


= K 




where K = y/—2mUQRo/Ih, and 

m 


V = 


V-S 


( 12 ) 


(13) 


K is typically of the order 1 — 5 for 1 > 3 (corresponding 
to the magic spin-orbit gaps). Since for the nucleon mass 
m « 940 MeV and V — S « 750 MeV: p = 1.25, it follows 
from Eq. (12) that for the nuclear system the ratio x is 
of the order 1 — 5, that is, in nuclei the energy splitting 
between spin-orbit partner states is comparable in magni¬ 
tude to the spacings between major oscillator shells. This 
is because of the near equality of the mass m and the po¬ 
tential V—S in nuclei. It should be noted than in the case 
of the pseudo-spin symmetry [T5] (V=-S), Eq. (131 yields 
7y=m/2V. An aspect that can be generalised to different 
quantum systems is the specific functional dependence of 
X on the ratio of the particle mass m and the effective 
potential whose gradient determines the spin-orbit force. 


III. OTHER SYSTEMS OF BOUND FERMIONS 


In the case of atomic systems the binding of an electron 
is determined by the Coulomb potential V{r) = —Za/r, 
where Z is the charge of the nucleus and the fine- 
structure constant a = 1/137. One can again perform 
the non-relativistic reduction of the Dirac equation for 
the electron (Eqs. and the resulting spin-orbit 

potential reads: 


yLS ^ 1 I dVjr) 

2M^ir)r dr 


( 14 ) 
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with 2A4{r) = 2m — V{r). In this case, of course, U{r) 
in Eq. ([^ is just the Coulomb potential. The energy 
spacing between successive levels with different principal 
quantum number n is proportional to whereas the 
first-order spin-orbit splitting is ~ a^. The ratio between 
the principal energy spacings and the spin-orbit splittings 
(fine structure) is much larger than in the nuclear case, 
that is ~ 1/a^ ~ 10 ^, known from the early seminal 
work of Sommerfeld M- The interesting fact is that, 
starting from Eq. (141, this ratio can again be expressed 
with the same functional dependence on 77 as in Eq. (12 1 . 
77 = m/V is now negative and, with m = 0.5 MeV and 
y(ro) = —2.72 X 10“^ MeV for the hydrogen atom and 
the Bohr’s radius tq, it follows that in the atomic case the 
characteristic value is 77 ^ —1/a^ « —2 x 10"^. For large 
absolute values of 77 , the expression Eq. (12) reduces to 
X ~ 1 /a^ « 10 ^, in agreement with the empirical value 
quoted above. The fine structure of atomic spectra thus 
becomes a limit of the spin-orbit rule (12|. It should be 
noted the validity of the spin-orbit rule in Coulomb-like 
systems is due to the 1/r behavior of the potential. The 
spin-orbit rules also applies to the case of ions having Z 
protons: in these systems the fine structure is known to 
behave as In the present approach the 7? factor 

comes from the Z one of Eq. (13) and the Z one from the 
rg/Z typical size of the ion. 

Figure displays the quantity representing the spin- 
orbit rule: 


X ~ 




(15) 


as a function of the ratio 77 between the mass of the par¬ 
ticle and the effective potential that determines the spin- 
orbit force in a given quantum system. As shown above, 
for nuclei 77 is slightly larger than one and, depending on 
the specific orbital, x lies in the interval 1 — 5. In atoms 
77 is negative and of the order of 10 ^ and, therefore, the 
characteristic value of a: is ~ l/a^ « lO'^. 

The main results of the spin-orbit rule quantatively 
apply to nuclei and atoms. It may be relevant to test it 
in a few other quantum systems although a more qual¬ 
itative agreement is expected, due to additional effects 
which would not be included in the present approach. 
For instance, experimental evidence indicates that in A 
hypernuclei the A-nuclear spin-orbit interaction is very 
weak compared to the strong spin-orbit interaction in 
ordinary nuclei [501 HI] ■ Available data show that the en¬ 
ergy spacings between A spin-orbit partner states are of 
the order of « 100 keV, although some microscopic mod¬ 
els also predict 2 MeV values in some cases [55]. Even 
though a quantitative study of the smallness of the A 
spin-orbit interaction necessitates a rather involved anal¬ 
ysis based on in-medium chiral SU(3) dynamics [53], one 
can qualitatively understand the reduction of the spin- 
orbit splitting in A hypernuclei already at the relativis¬ 
tic mean-field level [MlES]. A finite-density QCD sum- 
rule analysis nmi] indicates that, when compared the 
strong scalar and vector mean fields experienced by nu- 



FIG. 1: The ratio between the principal energy spacings 
and the spin-orbit splittings (fine structure) Eq. (15), as 
a function of the ratio 77 between the mass of the 
particle and the effective potential that determines the 
spin-orbit force in a given quantum system. 


cleons in ordinary nuclei, the corresponding self-energies 
of a A hyperon are reduced by a factor « 0.4 — 0.5, 
and even smaller if corrections from in-medium conden¬ 
sates of higher dimensions are taken into account. With 
TTTA = 1.12 GeV, this means that 77 in Eq. (13) is of the 
order 3 — 4, leading to larger x value than in nuclei, in 
qualitative agreement with empirical values of x in such 
systems. Since the energy spacing between major oscilla¬ 
tor shells in A hypernuclei (for instance, ^^G, a® 0, A°Ca), 
is of the order of « 10 MeV [20], the empirical ratio x 
ranges from a few units to a few dozens. Our mean-field 
estimate is closer to the small values of x but, as shown 
in Ref. |5S], a considerable additional reduction of the 
splitting between A spin-orbit partner states arises as the 
effect of the hypernuclear tensor coupling. As explained 
above, a quantitative explanation of the small spin-orbit 
splitting in hypernuclei must include beyond-mean-field 
effects. 

An interesting example of quantum systems that are 
governed by the strong interaction but exhibit negative 
values of 77 , are quarkonia such as, for instance, char- 
monium cc and bottomonium bb. The center-of-mass is 
different, compared to the present one-body approach, 
but this only generates a global scaling factor of 2 on the 
energy positions, due to the reduced mass. Since we are 
dealing with order of magnitudes, and most importantly 
with energies ratios, it is relevant to check if the present 
approach applies in the case of quarkonia. The use of 
potential models for these system can be justified by the 
fact that the bottom and charm masses are large in com¬ 
parison to the typical hadronic scale of QGD. Most phe¬ 
nomenological approaches to the dynamics of two heavy 
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quarks interacting through a potential are variants of the 
Cornell model [26], which consists of a superposition of 
the one-gluon-exchange that leads to a Coulomb-like at¬ 
tractive vector potential at short distances: 


^ 4 a^(r) 

3 r 


(16) 


plus a scalar linear confining potential S = or, with 
a « 0.18 GeV^. Therefore, 1^ — 5” is negative in the 
quarkonia case. The masses of the c quark b quark are 
rric « 1.27 GeV and rub « 4.2 GeV, respectively. At ra¬ 
dial distances that correspond to the mean-square radii 
of the quarkonia states the Coulomb-like attractive vec¬ 
tor potential Eq. (16 1 , with a depth of the order of —1 
GeV, dominates over the scalar potential [23 . The mass 
of the lowest charmonium state is « 3 GeV, while that 
of the lowest bottomonium state is « 9.5 GeV. The en¬ 
ergy spacing between IS' and 2S states is of the order 
of 600 MeV for charmonia, and 560 MeV for bottomonia 
[2E1I1S]- The fine splittings between IPj states are less 
than 100 MeV for the charmonium (AM 21 = 45.6 ± 0.2 
MeV and AMio = 95.3 ± 0.4 MeV), and less than 40 
MeV for the bottomonium (AM 21 = 19.4 ±0.4 MeV and 
AMio = 33.3±0.5 MeV). This means that, in the case of 
quarkonia, the empirical ratio between the energy spac- 
ings characterized by the principal quantum numbers and 
the fine spin-orbit splittings is of the order of 5 to 10, in 
qualitative agreement with Fig. 1. 

The particular dependence of this ratio on rj, shown in 
Fig. HI displays interesting features. For positive values 
of rj, states for which the orbital angular momentum and 
spin are aligned are found at lower energy with respect 
to states for which the orbital angular momentum and 
spin are anti-aligned (nuclei, hypernuclei), whereas the 
opposite energy ordering is found for negative rj (atoms, 
quarkonia). The ratio x (Eq. (HU)) diverges at 77 = 0 , 
that is, in the limit of massless particle. 

An interesting situation is found in the vicinity of 
77 = 1 / 2 , for which x = 0. More precisely, 77 values close 
to 1/2 are in the validity domain of the approximation 
performed to derive the present spin-orbit rule. This oc¬ 
curs when the mass of the particle is close to {V — S)/2, 
and the energy spacings between states are characterized 
by very large spin-orbit coupling (giant LS). Such states 
could be obtained in particular cases for which one would 
be able to choose the strength of the effective potential 
whose gradient determines the spin-orbit force. In a hy¬ 
pothetical case this limit is approached when the effec¬ 
tive potential becomes very deep as it would occur, for 
instance, for bound antibaryon-nuclear systems which, 
in addition to ordinary nucleons, contain antibaryons 
{B = p,A,...) (cf. Figs. 10 and 21 of Ref. [30]). 
Two cases can be considered. For the antibaryon spec¬ 
trum the spin-orbit splitting remains small, of about few 
hundreds keV m- The corresponding V-S potential is 
smaller than in the nucleon case because V has the op¬ 
posite sign in the case of antibaryons |31| and, therefore, 
the value of 77 becomes larger (Fig. H])- This is both the 


case for the antiproton and anti-A nuclei [221132] : in the 
last case a good quantitative agreement is found for the 
relative amplitude of the spin-orbit splitting predicted by 
the spin-orbit rule (x^ 40). 

The case of nucleons in antibaryon-nuclei is closer to 
the giant LS state: the gradient of the potential obtained 
from self-consistent calculations increases and the spin- 
orbit splitting not only gets larger, but displays a pro¬ 
nounced increase with respect to the spacing between 
major shells (Fig. [^. However, such systems have not 
been observed and this is just an illustrative example, 
with antibaryon-nucleus potentials subject to large un¬ 
certainties [30]. A similar effect, although much weaker, 
is predicted for single-nucleon spectra in A-hypernuclei 
[33] • Even though the presence of a A induces only a frac¬ 
tional change in the central mean-field potential, through 
a purely relativistic effect it increases the spin-orbit term 
in the surface region, providing additional binding for 
the outermost neutrons. To emphasize this region of gi¬ 
ant spin-orbit splittings, in Fig. [l] we also plot the inverse 
of X. For negative values of 77 the ratio x does not vanish, 
but again one finds a region of small absolute values of 77 
for which the splitting between spin-orbit partner states 
is comparable in magnitude to the energy spacing be¬ 
tween successive levels with different principal quantum 
number. 


IV. SUMMARY 

In conclusion, we have investigated an interesting effect 
of spin-orbit coupling in systems of bound fermions (elec¬ 
trons, quarks, nucleons, hyperons). By performing the 
usual non-relativistic reduction of the Dirac equation, the 
single-particle mean-field equation takes a Schrddinger- 
like form which, in addition to the confining potential, 
exhibits the spin-orbit potential explicitly. Starting from 
the relativistic mean-field approximation in the nuclear 
case, we have derived an analytic expression for the ratio 
between the energy spacings characterized by principal 
single-particle quantum numbers and the energy split¬ 
ting of spin-orbit partner states. This quantity explic¬ 
itly depends on the ratio of the particle mass and the 
effective potential whose gradient determines the spin- 
orbit force (cf. Eq. ([T^). In nuclei the near equality 
of the mass of the nucleon and the difference between 
the large repulsive vector and attractive scalar potentials 
explains the fact that the spin-orbit splittings are com¬ 
parable to the energy spacing between major oscillator 
shells. The same universal functional form also applies 
to other bound quantum systems, regardless of whether 
the spin-orbit potential originates from the strong or elec¬ 
tromagnetic interaction, and for which the ratio between 
the principal energy spacings and the spin-orbit split¬ 
tings is orders of magnitude larger than in the nuclear 
case. For a given particle mass this ratio could be al¬ 
tered, in principle, by modifying the potential that gen¬ 
erates the spin-orbit coupling effect. When this ratio 
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is rj « 1/2, our study predicts the occurrence of giant 
spin-orbit splittings, that is, a single-particle excitation 
spectrum dominated by large energy spacings between 
spin-orbit partner states. 

The extension of the derivation of the spin-orbit rule 
to exotic nuclei, where the densities are more diffuse shall 
be undertaken in a forthcoming work. This includes halo 
nuclei involving coupling to the continuum [34H37| . It 
would also be relevant to derive a pseudo-spin rule, pro¬ 


viding the typical magnitude of the related degeneracy 
raising due to the breaking of the pseudo-spin symmetry 
[HUSH], with respect to the main shell structure. 
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